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I. INTRODUCTION 

Catastrophes jij are at the heart of many fascinating 
optical phenomena. The most prominent example of such 
a catastrophe is the rainbow. Light rays from the Sun 
enter water droplets floating in the air. After two refrac- 
tions and one reflection inside each drop the rays reach an 
observer. Above a critical observation angle no rays ar- 
rive, whereas below the angle two rays strike the observer. 
A bright bow, the rainbow, appears at the critical angle, 
because here the cross section of light rays diverges [||. 
(The critical angle depends on the refractive index that 
varies with the frequency of light in dispersive media such 
as water, giving rise to the rainbow colors.) The direction 
of a light ray is proportional to the gradient of the phase. 
The rainbow thus represents a singularity of a gradient 
map, a catastrophe in the sense of Thom |^| and Arnol'd 
j| . Structurally stable singularities of gradient maps fall 
into distinct classes, depending on the number of control 
parameters involved (3|i|]. Structural stability is the key 
to Nature's way of focusing light [^) in the caustics cre- 
ated by ray catastrophes. Yet the wave nature of light 
smooths the harsh singularities of rays. Simultaneously, 
characteristic interference effects appear. For example, 
the pairs of light rays below the rainbow create a delicate 
pattern of supernumerary arcs jl[] that are visible under 
favorable weather conditions (when the floating droplets 
are nearly uniform in size Q). Every class of diffraction 
catastrophes generates its distinct interference structure 

Catastrophe optics describes the wave properties of ray 
singularities. In the hierarchy of physical concepts, wave 
optics refines and embraces ray optics, and quantum op- 
tics rules above wave optics. So, what would be the 
quantum effects of wave catastrophes (6|? First, what 
are quantum catastrophes? It might be a good idea to 
begin with an example, the black hole (7j. When a star 
collapses to a black hole an event horizon is formed, cut- 
ting space into two disconnected regions. Seen from an 
outside observer, time stands still at the horizon, freez- 
ing all motion. A light wave would freeze as well, propa- 
gating with ever-shrinking wavelength. In mathematical 



terms ||, a monochromatic light wave of frequency lo 
oscillates as Q(r — r s )(r — r s ) lM when the radius r ap- 
proaches the horizon r s , with fi — 2r s lo/c. A logarith- 
mic phase singularity will develop. Potential quantum 
effects of such a wave singularity are effects of the quan- 
tum vacuum. The gravitational collapse JtJ of the star 
into the black hole has swept along the vacuum. The vac- 
uum thus shares the fate of an inward-falling observer. 
Yet such an observer would not notice anything unusual 
at the event horizon. In mathematical terms, the vac- 
uum modes are analytic across the horizon |8|,^j . On the 
other hand, the modes perceived by an outside observer 
are essentially non-analytic, because they vanish beyond 
the horizon where the observer has no access to. Con- 
sequently, the observer does not see the electromagnetic 
field in the vacuum state. Instead, the observer notices 
the quanta of Hawking radiation JTo| with the Planck 
spectrum (e 27rM — The quantum vacuum does not 

assume catastrophic waves of the type 0(r — r s )(r — r s )* M , 
hence resolving so the associated wave singularity and, si- 
multaneously, generating quantum radiation with a char- 
acteristic spectrum. At the heart of such a catastrophe 
lies a time-dependent process, for example the gravita- 
tional collapse in the case of the black hole §. The 
process has disconnected the spatial regions where waves 
can propagate and has created a logarithmic phase singu- 
larity at the interface. Any time-dependent phenomenon 
will generate some radiation, as long as the process lasts. 
In remarkable contrast, a quantum catastrophe creates 
quanta continuously. 

This paper describes the theory behind a recent idea 
[11 1 to generate a quantum catastrophe of slow light 
[12 p^| , p^7| -[2l|] . An experiment is proposed based on 
Electromagnetically-Induced Transparency (EIT) Jl7| . 
In EIT a control beam determines the optical properties 
of slow light in a suitable medium. Changing the inten- 
sity of the control light from a uniform to a parabolic 
profile creates s slow-light catastrophe Oj, see Fig. 1. 
This catastrophe resembles the event horizon of a black 
hole but also shows some characteristic differences. In 
Section II we put forward a rather general phenomeno- 
logical quantum field theory of slow light. Appendix A 
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justifies the theory for the most common method to slow 
down light using EIT JlTj . In Section III we address the 
specific theory of the slow-light catastrophe |11|. Ap- 
pendix B contains some estimations that are relevant to 
the experimental aspects involved. Section IV summa- 
rizes the results and draws a further vision of quantum 
catastrophes. 




FIG. 1. Schematic diagram of the proposed experiment. A 
beam of control light with intensity I c generates Electromag- 
netically-Induced Transparency [15] in a medium, strongly 
modifying its optical properties for a second field of slow light. 
When an initially uniform control intensity is turned into the 
parabolic profile shown in the figure, the slow-light field suf- 
fers a quantum catastrophe. To slow-light waves, the interface 
Z of zero control intensity cuts space into two disconnected 
regions and creates a logarithmic phase singularity, in analogy 
to the effect [7] of an event horizon [6] . The quantum vacuum 
of slow light cannot assume such catastrophic waves. In turn, 
pairs of slow-light quanta, propagating in opposite directions 
away from Z, are emitted with a characteristic spectrum. The 
waves shown below the intensity profile refer to the emitted 
light with the modes vjr and wl of Eq. (fed). 



II. THE MODEL 

Electromagnetically-Induced Transparency (EIT) |l7]] 
has served as a method to slow down light significantly 
[ p2[ or, ultimately, to freeze light completely [jl3] [Tg[ . 
Like other successful techniques, EIT is based on a simple 
idea |l7j] : A control beam of laser light couples the upper 
levels of an atom, and, in this way, the beam strongly 



modifies the optical properties of the atom. In partic- 
ular, the coupling of the excited states affects the tran- 
sition from the atomic ground state to one of the up- 
per states, i.e., the ability of the atom to absorb probe 
photons with matching transition frequency. Destructive 
quantum interference between the paths of the transi- 
tion process turns out to eliminate absorption at exact 
resonance ]l7| ]. A medium composed of such optically- 
manipulated atoms is transparent at a spectral line where 
it would otherwise be completely opaque. In the vicinity 
of the transparency frequency the medium is highly 
dispersive, i.e. the refractive index changes within a nar- 
row frequency interval. In turn, probe light pulses with 
a carrier frequency at luq travel with a very low group 
velocity v g [|l8| . The intensity I c of the control beam de- 
termines the group velocity of the probe pulse, as long 
as I c is stronger than the probe. Paradoxically, the lower 
I c is the slower moves the pulse, which, however, is only 
possible when the electronic states of the atoms follow dy- 
namically the control field fl9fl, causing the probe-light 
intensity to fall accordingly |l9| . In this regime light 
freezes when I c vanishes 1L| fl(j| . 

The theory of EIT |l7]] often employs an atomic three- 
level scheme: Two levels account for the excited states 
coupled by the control field and one level represents the 
ground state, see Appendix A. In reality, atoms are more 
complicated and, when details matter, an accurate de- 
scription involves the full atomic sub- level structure [^0) . 
Here we put forward a phenomenological quantum theory 
of slow light that is rather independent of the microscopic 
mechanisms used in practice. We assume only that the 
slow-light medium is transparent with a real susceptibil- 
ity that depends linearly on the detuning from luq. In 
EIT our model is restricted to the narrow transparency 
window around luq that can maximally reach the natural 
line width of the atomic transition. Appendix A shows 
that our model agrees with the three-level theory of EIT. 
Our theory is simple enough to be treated analytically 
and yet sufficiently complex to capture the essence of 
slow-light quanta. We postulate an effective Lagrangian 
J*? , show that Jz? is consistent with the known dynamics 
of EIT within the validity range of our model, calculate 
the energy and quantize the field of slow- light polaritons. 
Our approach has the additional advantage that it may 
be applicable to other mechanisms plf for creating slow 
light that do not rely on EIT. 



A. Lagrangian 

We characterize slow light by a real scalar field ip ignor- 
ing the polarization. The optical field (p shall be given in 
units of the vacuum noise p2[ . For simplicity, we assume 
uniformity in two spatial directions, x and y, and regard 
the optical field as a function of time t and position z. 
Throughout this paper we denote partial time derivatives 
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by dots, by differential operators dt or simply by suf- 
fixes t, whatever happens to be most convenient. Spatial 
derivatives are denoted by dashes, differential operators 
d z or suffixes z. We postulate the effective Lagrangian 
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We see in the next subsection that the real parameter 
function a determines the group velocity, and hence we 
call a group index. In EIT jlTj the parameter a is in- 
versely proportional to the control-field intensity, 



a(t, z) 



I c (t, z) 



(2) 



with a coupling strength n that is proportional to the 
modulus squared of the atomic dipole-transition matrix 
element and to the number of atoms per unit volume. 
Without the EIT medium present, Jz? is the Lagrangian 
of a free electromagnetic field, ^(E 2 — c 2 B 2 ) with fixed 
polarization. We see that ip is related to the electric field 
strength E in SI units by 
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where to denotes the frequency of light and £o is the elec- 
tric permittivity of the vacuum. 

In order to motivate the Lagrangian ([!]) we consider 
the corresponding Euler-Lagrange equation p3| , 
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that leads to the wave equation 
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This is the propagation equation of slow light based on 
the traditional three-level model, see Appendix A, but 
the equation holds on more general grounds: Equation 
(||) describes light in media with linear spectral suscepti- 
bility. Assuming that the optical field oscillates at much 
shorter time and length scales than any variations of a, 
we replace idt by the frequency to and —id z by the wave 
number k. We arrive at the dispersion relation 
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which, in the positive transparency window near u>o, 
agrees with 







and the linear spectral susceptibility 
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At a later stage we need to consider negative frequencies 
as well. To verify that the dispersion relation (^) is also 
valid in the corresponding negative transparency window 
we utilize the general property of a spectral susceptibility 
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We see that as long as the frequency of the probe light 
lies within the transparency windows of EIT, the La- 
grangian (^) reproduces the typical linear slope of the 
spectral susceptibility. In fact, up to a trivial prefactor, 
_Sf is the only Lagrangian that is quadratic in the field 
and its derivatives and that is consistent with the spec- 
tral susceptibilities (||) and (|l0|) . Therefore we regard Jzf 
as a suitable Lagrangian for slow light. 



B. Dynamics 

The ability to freeze light by turning off the control 
field depends crucially on the dynamics of the process. 
Consider a time-dependent group index a without sig- 
nificant spatial variations. In this case slow light is dom- 
inated by oscillations within an optical wave length and 
an optical cycle. We express the wave as 
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with the slowly varying electric-field amplitude £ in SI 
units. We approximate 
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and get from the wave equation (^J) 



-2iw(l + a)£ 

w ((1 + a)oj 2 + iusa - <?k 2 + 2ikc 2 d z - aco 2 ) £ 
= (2iujcd z + iu>a + a(u> 2 — w )J £ 
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when the carrier frequency to is equal to the transparency 
resonance loq. We apply the relation (||) between the 
control-field intensity I c and the group index a, write I c 
as the square of the field strength £ Cl and obtain, finally, 
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£ + cE ' = -a£ - ^£ = d t ^r . 



(14) 



This is exactly the propagation equation of slow light 
in the adiabatic and perturbative limit (Eq. (9) of Ref. 
|jl9f with the Rabi frequency f2 being proportional to £ c ) . 
Consequently, the Lagrangian (|l|) has codified naturally 
the correct dynamic regime including the ^a£ term that 
describes reversible stimulated Raman scattering iQ . 

In order to understand the principal behavior of or- 
dinary slow-light pulses, consider the case of a spatially 
uniform yet time-dependent group index. Equation Jl4| ) 
has the solution 

£(t, z) = £ (z - Jvgdtj \jv g jc (15) 



in terms of the velocity yA 



1 + a 



(16) 



We see that the pulse envelope £ propagates with the 
speed v g called group velocity. When the group velocity 
changes in time, the intensity \£ | 2 reacts proportionally. 
The ratio between the control (^) and the pulse intensity 
(|l5|), (J c + k)/ \£q\ 2 , remains large, even in the limit of a 
vanishing control field when v g vanishes as well, as long as 
k is large (for a sufficiently dense medium) . The spectral 
spread of a pulse is reduced by v g /c and a standing pulse 
oscillates just with the carrier frequency. In this regime 
slow light does not leave the transparency window of EIT 
E3 . One can freeze light without losing control [p3UT3] . 



C. Energy 

After having gained confidence in our field-theoretical 
approach, we use the Lagrangian ([!]) to calculate the en- 
ergy balance of slow light. According to Noether's theo- 
rem [ p3| we obtain the energy density 

/=^-^=| ((1 + a)^ + cV 2 + ac V) (17) 



and the energy flux (Poynting vector) 

Jr = -z— (p = —he (pip . 
dip' 



(18) 



The energy balance It + P z is then, as a consequence of 
the wave equation (|B|), 



(19) 



Temporal changes in the control field, modifying the 
group index (^J), do not conserve energy. In fact, the ex- 
periment JT~3[ ] indicates that the control beam can amplify 
light stored in an EIT medium with zero group velocity. 



In the experiment |f|], slow light enters the EIT sam- 
ple and is then frozen inside by turning off the control 
field. Switching on the control releases the stored light. 
The pulse emerges with an intensity that depends on the 
control field and that may exceed the initial intensity, 
in agreement with Eq. (|l5|). Clearly, this phenomenon 
is only possible if energy is indeed transferred from the 
control beam to the probe light. 



D. Polar itons 

Finally, we realize the full potential of the Lagrangian 
(|l[) in setting up an effective quantum theory of slow 
light. The classical canonical momentum density of the 
field p is @J2§ 
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We quantize the field by regarding ip and SJC/Scp as Her- 
mitian operators ip and 7T, respectively, with the canoni- 
cal commutation relations l|23],|25|l 



[$(t,z),<p(t, Z')] = [7t{t,z),7T(t,z')} =0, 

[<p(t ) z) ) Tt(t,z')] = i%6(z-z'). (21) 

Let us decompose the field (p into modes with dimension- 
less mode indices q 



^(M) = / (a q u q (t, z) + a\u*(t, z)) dq . 



(22) 



In order to guarantee that <p satisfies the wave equation 
(||) the mode functions u q are required to obey Eq. (|^) 
as well. The u q shall be normalized according to 

(u q ,u q/ ) = 6(q- q') , (u q ,u*,)=0, (23) 

with the Klein-Gordon-type scalar product (2^j2j| 

((pi,(p 2 )=i ((pi(p2-(p*ip2j(l + (x)dz. (24) 

J — 00 

The scalar product is chosen such that it remains con- 
stant during the propagation of (pi and tp2, 



dt (¥>i,P2) 
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Using these postulates and definitions we calculate the 
commutation relation of the mode operators 
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In a similar way we prove that 



0. 



(27) 



Consequently, and in agreement with the spin-statistics 
theorem [^3|, slow light consists of bosons. Let us ex- 
press the total energy ([l8]) in terms of the annihilation 
and creation operators a q and d q . Consider the case of 
a stationary group index a when the total energy is con- 
served. We obtain after partial integration, via the wave 
equation (|) for the field operators, 



f 0Q dz 
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((1 + a){d t <p) 2 - c 2 ip d 2 z + a ujI<p 2 ) dz 
({d t cp) 2 -tpdfip) (l + a)dz. (28) 



We employ the mode expansion ( p2|) w ith the norm ( |23| ) 
with respect to the scalar product (|24|), use the commu- 
tation relation (E6[), and find, finally, 



T w dz = I hcu [ a\a q + I ) 



(29) 



Consequently, the annihilation and creation operators a q 
and a q refer indeed to the energy quanta of slow light. 
The quasiparticles of light in a dielectric medium are 
called polaritons in analogy to the optical excitations in 
a solid |20| . They combine characteristic features of free 
photons with the properties of the dipole oscillations of 
the atoms constituting the medium. When light is slowed 
down photons turn into atomic excitations that, after ac- 
celeration, may emerge as photons again |l^,[l4|]. The 
polariton picture contains implicitly the correct book- 
keeping of the photonic and atomic features of light in 
linear media. 



One can use similar arguments as in the subsection on 
the dynamics of slow light to prove that the proposed po- 
lariton theory is consistent with the adiabatic three-level 
model [^9) . Yet our approach is not restricted to a regime 
dominated by spatial oscillations of the form cxp(io;z/c), 
see Eq. (^Tj), which is in essence the regime of geomet- 
rical optics |27). One can easily relax this unnecessary 
constraint of the three-level theory, because the electro- 
magnetic response of an atom is local, as long as the wave 
length of light is large compared with the size of the atom. 



The adiabatic polariton theory |19| is, as ours, restricted 
to a narrow frequency range with respect to time, i.e., 
to the transparency window of the electromagnetically- 
manipulated medium. Yet this restriction in frequency 
does not exclude rapid spatial oscillations beyond the 
scale of the wave length in vacuum, 2ttc/uj. We will see 
in the next section that such oscillations occur near a 
spatial singularity of the group index. In this situation 
our quantum field theory of slow-light polaritons turns 
into a perfect tool for analyzing the quantum physics of 
a wave catastrophe. 



III. THE CATASTROPHE 

Imagine that the control beam illuminates the EIT 
medium from above. Initially, the control intensity is uni- 
form, but then the control light develops a dark node that 
continues as an interface Z of zero intensity I c through a 
part of the medium. One could use computer-generated 
holograms to achieve this situation, similar to the gener- 
ation of Laguerre-Gaussian beams pgj ]. Suppose that the 
interface is sufficiently flat and cuts deep enough into the 
medium to justify our model. (We have assumed unifor- 
mity in the two spatial directions x and y that are par- 
allel to the interface. Most probably, uniformity over a 
few wave lengths would suffice.) Near a zero in the inten- 
sity, the control field strength must grow linearly. Con- 
sequently, I c depends quadratically on z and, according 
to the relation (^), the group index a forms a quadratic 
singularity where the group velocity (|l6| ) vanishes, 



(30) 



The parameter a sets the scale of the group-index pro- 
file. We assume that the spatial profile (|3FJ) of the group 
index extends over a sufficiently long range. For simplic- 
ity, we consider a one-dimensional model where the slow 
light propagates in z direction only. Appendix B general- 
izes our results to a realistic three-dimensional situation. 
Apart from these idealizations and from the physics cap- 
tured in our Lagrangian ([!]) , we make no further approx- 
imations. 
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A. Horizon 



v = m 



(37) 



At a node of the control field the group velocity of the 
probe light vanishes. Therefore we would expect that 
no wave packet of slow light can pass the interface Z. 
Consider slow light subject to the wave equation (§) with 
the group- index profile (30). We represent a wave ip(t, z) 



<p = ^/z(f> 



(31) 



and find 



[z 2 d 2 t + a 2 (d 2 + u; 2 ) - c 2 (zd z zd z - i)] cf> = 0. (32) 

We could multiply 4> by step functions 0(±z) and still 
get solutions of the wave equation ([32]), because 

zd z Q(±z)</>(z) = e{±z) zd z <j)(z) ± <j){z) z 5{z) 

= e{±z)zd z <j)(z). (33) 

Consequently, waves on different sides of Z are indepen- 
dent from each other. As long as slow light is concerned, 
the interface Z has cut space into two disconnected parts. 

Consider monochromatic probe light oscillating with 
frequency u. In this case the wave equation ( |32| ) reduces 
to Bessel's differential equation p9| with the index 



a 2 (fc 2 - k 2 ) 



(34) 



Here k abbreviates u>/c and ko refers to loq/c. So, in 
mathematical terms, the monochromatic waves of slow 
light are products of a square root with Bessel functions 



ip = \fz J± u (kz) e 



-iut 



(35) 



Depending on the detuning of the frequency lo with re- 
spect to the exact transparency resonance u>q, two cases 
emerge. First, when 4a 2 (k 2 — k 2 ) is below unity the Bessel 
index v is real. We apply the asymptotics of the Bessel 
functions for large and positive arguments p [p9[, 
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We see that in the far field the light waves with real 
Bessel indices are in a perfectly balanced superposition 
of incident and emerging plane waves. In other words, the 
light is totally reflected away from the interface of zero 
group velocity, similar to the reflection of radio waves at 
the Earth's ionosphere [3C| ]. 

A more interesting scenario appears in the second case 
when the light is sufficiently blue-detuned to evoke an 
imaginary Bessel index 



This regime can be reached by adjusting the gradient of 
the control field that plays a decisive role in Eqs. (||) and 
( pp| ) . The smaller the gradient is the larger is a 2 . There- 
fore, a sufficiently small control-field gradient gives rise 
to an imaginary Bessel index v, even within the narrow 
transparency window in frequency lo = kc. On the other 
hand, as we will see at a later stage, the gradient should 
be as large as possible for producing a maximal quantum 
effect. When the Bessel index v is imaginary the inci- 
dent and the reflected waves are not balanced anymore, 
as we see from the asymptotics (^6|) of the Bessel func- 
tions. Only a fraction of the incident wave is reflected 
and the rest must be transmitted somewhere. To find 
the transmitted wave, we focus on the behavior of the 
Bessel functions for small arguments. We use the first 
term in the power series |29(| , 
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£ = kz, C . 



(38) 



Therefore, near the interface of zero group velocity the 
light waves are proportional to 



^+1/2 = ^ exp (l/ilnC) = ^ exp( . 5) 



(39) 



The logarithmic phase S reduces dramatically the wave 
length near z = 0, because 



2tt 2tt 
A = — = — z. 
b z fi 



(40) 



The transmitted wave thus freezes in front of the interface 
of zero group velocity. 

We regard a process that creates an interface where 
waves separate and develop a logarithmic phase singu- 
larity as a wave catastrophe. The gravitational collapse 
of a star into a black hole has created an horizon that 
cuts space into two disconnected parts as well [Q. Close 
to the event horizon, an outside observer would see a 
similar behavior of waves || . All motion freezes near the 
horizon of the hole which, therefore, has also been termed 
frozen star . In view of this analogy we regard the in- 
terface of zero group velocity Z as the optical analog of 
a horizon. 



B. Modes 

Consider a superposition of sufficiently blue-detuned 
slow- light waves ([35]) with imaginary Bessel indices (|37]) . 
To analyze the quantum effects of the horizon on polari- 
tons, we must turn the waves into modes i.e., we 

must normalize the wave functions according to the scalar 
product (p4|). Waves (^) with different wave numbers 
k are orthogonal to each other and possess a continu- 
ous spectrum. Hence they should be normalized to delta 



G 



functions. We employ the ratio fe/feo as the dimensionless 
mode index q that occurs in the scalar product ( p4| ) . The 
normalization factor is entirely determined by the way in 
which the norm (^4) diverges to reach the S singularity 
p2[ . Consequently [32j, we can ignore all finite, converg- 
ing contributions to the normalization integral given by 
Eqs. (pj), m, and ©, 



(91,92) =(wi+w 2 ) / Jh^kxz) J ±l ^ 2 {k 2 z) . 
Jo 



x [ 1 + — ) dz . 
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The integral (|4l]) diverges both at z = and at z = 
00. We account for the two divergences seperately in the 
integrals Iq and Jqo, 



(91,92) = 2w (Iq + Ioc) 



(42) 



Ignoring any convergent contributions to Iq, we cut off 
the integral at some small zo, regard z(l + a 2 / z 2 ) as a 2 /z, 
and use the asymptotics ( |38| ) near the origin, 
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sinh(7r/z) — S(fj,2 — fJ-i) 
fc _1 sinh(7r/x) — fc 2 ) . 



(43) 



In the last step we have utilized that fj,(d/j,/dk) equals 
a 2 k. Furthermore, we have used Eq. 1.2.(6) of Ref. |3^] 
for the gamma function a;! = T(x + 1). Let us address 
the other integral, 
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2Vfc 



J± itll (k x z) J±i M2 (k 2 z) z dz 
exp y— ik\z ± fii — +i 
+ exp {+ik\Z TfJtiZ ~ 



exp 



-tkoZ ± fro 1 — 

p 2 4 



/ ., 7T ,7T 

+ exp \-ik 2 z =F ^2^- + «j 
exp (i(fc 2 - h)z ± (^1 + ^2)7^) 

/ 7T 

+ cxp(^i(fci -A: 2 )2:=F (Mi + ^2)2 



G?2 



1 

2?rfc 



= k 1 cosh(7r/i) <5(fci — k 2 ) ■ 



dz 
(44) 



Combining the two integrals (^) and (^4|) gives in total 



(^,^2) =2ce«»5{k 1 



(45) 



Consequently, in the spatial region right from the horizon 
(z > 0) the normalized wave functions are 

u% = e -T/./2 ^ j (fcz) e"*-* . (46) 
V 2c 



Left from the horizon (z < 0) we chose for convenience, 

ut(z)=ul(-z). (47) 

The step function in the definition ( f46| ) guarantees 
that the u modes on different sides of the horizon do not 
overlap and, therefore, they are automatically orthogonal 
to each other. However, the ± degenerated waves on the 
same side are not orthogonal. In fact, we obtain along 
similar lines as in the normalization procedure, 



"R, "-RJ 



%1 - 92) 



(48) 



Yet we construct easily the orthogonal partners w to 
the u T waves, 



1 



w 



R 



l Rl ' 



wl(z)=wl(-z). (49) 
We chose as the orthonormal basis in the mode expansion 

(H) 

w R = w^, u R = u R , Wl=w2, u L =u~l . (50) 

Finally, to find an interpretation of the w modes, we use 
the asymptotics ([56]) of the Bessel functions and get for 
\z\ — > 00 on the appropriate sides of the horizon, 



WR 



1 - e 



-2iTfi \ I/ 2 



1 — e~ 27rM 
47TC uj/ujo 



^ exp ( 



LO IT 

+i— z — iujt — i — 
c 4 



1/2 



[-1-, 
V c 



exp ( — i— z — iu>t — i— ) . (51) 



The asymptotics ([51]) shows that the w modes turn into 
plane waves propagating away from the horizon. In other 
words, the w modes are the ones that reach an external 
photon detector. 



C. Analyticity 

The stationary modes ( |50| ) of catastrophic slow light 
are severely non-analytic. The modes vanish on either 
the left or the right side of the horizon with a charac- 
teristic essential singularity as a precursor. Waves near 
the event horizon of a black hole suffer a similar fate [|| . 
Seen from an outside observer, the waves freeze near the 
Schwarzschild radius r s with an essential singularity of 
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the type (r — r s ) v where fi — 2r s ui/c ||. Yet an ob- 
server falling into the hole would see little difference in 
waves near the horizon and could pass the point of no re- 
turn without noticing. Like the inward-falling observer, 
the quantum vacuum flows towards the central singular- 
ity of the hole and, similarly, the horizon should not be 
a special place for the vacuum either. In mathematical 
terms, the wave function of the vacuum is analytic |^||5). 
Consequently, the modes seen by the outside observer 
must not be in their vacuum states. In fact, they carry 
the quanta of Hawking radiation |M . The history of the 
hole formation during a gravitational collapse turns out 
to be responsible for the analyticity of the vacuum |^],[| . 
Inspired by the analogy between a black hole and our 
slow-light catastrophe, let as consider the history of our 
horizon. 

Suppose that the group index a was initially a largely 
uniform ao- Then, by tuning the control field, the group 
index develops a quadratic singularity, for example as 
a = a 2 (t)/(z 2 + b 2 (t)) with a 2 (-oo) -> oo, 6 2 (-oo) -> oo, 
a 2 (-oo)/6 2 (-oo) -> a , and finally b 2 {+oo) -> 0. The 
details of the process do not matter. Immediately af- 
ter a time t c the group index will possess the quadratic 
singularity that we are studying, creating a slow-light 
catastrophe. Given a uniform group index as the initial 
condition and no slow-light injected, the polariton vac- 
uum occupies initially packets of plane waves that we can 
sort into right- or left-traveling waves, 

ip± — J A±(uj) exp (^ z ^~ \j 102 + a o( U! ' 2 ~ w o) — iut^j 
xdoj. (52) 

Regarded as a function of complex z, the ip± waves are 
analytic in the upper (+) or lower (— ) half plane, respec- 
tively, because here the integral ( |52| ) converges. When 
the control field creates a horizon the vacuum modes 
must follow the wave equation (|j|). We assume that a 
is analytic apart from poles. Consider closed contour in- 
tegrals in either one of the half planes. We obtain from 
the wave equation (^), 

dt <j) (1 + a)tp t dz = c 2 <j> tp zz dz — luq <j) aip dz . (53) 



Due to the analyticity of the initial wave packets (152 
§ ip dz and § tp t dz were initially zero. Equation ( |53| ) in- 
dicates that both integrals remain zero, as long as a is 
analytic. At single poles of a we get ipt = — w 2 ip, which 
cannot generate a singularity. Higher poles of a do not 
contribute to the closed contour integrals. Consequently, 
the vacuum wave functions are always analytic in z. 



\ z 









FIG. 2. Space-time diagram of a slow-light catastrophe. 
The figure illustrates the fate of a wave packet ip(t, z) that ex- 
periences the formation of the horizon Z. Initially, the packet 
oscillates with positive frequencies in time t and propagates 
from the left to the right in space z. The horizon cannot gen- 
erate negative frequencies in the reflected light, apart from a 
brief burst that we neglect. On the left side of Z we thus re- 
gard ip(t, z) as analytic in t on the lower half of the complex t 
plane. Furthermore, <p(t, z) is analytic in z on the upper half 
plane throughout the history of the wave packet, because the 
process (jE|) conserves analyticity. Yet ip(t, z) is not analytic 
in t on the other side of the horizon, as the solution (|H^) indi- 
cates. Here waves with negative frequencies are continuously 
peeling away from the horizon, corresponding to a stationary 
creation of slow-light quanta. 

Consider the analytic properties of the vacuum waves 
with respect to time. Picture a wave with positive fre- 
quencies incident from the left, see Fig. 2. After the 
catastrophe a part of the wave may freeze at the horizon 
and the rest is reflected. If the wave happens to arrive 
during the formation of the horizon a brief burst of light 
with negative frequencies may be generated. However, 
in the stationary regime we are interested in, the re- 
flected light contains always positive frequencies. There- 
fore we regard the wave function ip(t, z) on the left side of 
the horizon as analytic in t on the lower complex plane. 
Given an analytic signal ip(€) at some z we can propa- 
gate it in space according to the wave equation (|^) with- 
out loss in analyticity, but we cannot pass the horizon, 
because here a is singular. It is therefore conceivable 
that beyond the horizon ip(t) is not analytic anymore. 
In other words, tp(t) may contain negative frequencies. 
In fact, we show in the next subsection that negative 
frequencies in time are unavoidable for not running into 
conflict with the analyticity in the spatial coordinate z. 
Waves propagating to the right are analytic in z on the 
upper half plane and they must have originated from a 
wave incident from the left. The analytic properties of 
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the vacuum waves in space are thus connected with their 
analytic properties in time. Analyticity in the upper half 
of the z plane is linked to analyticity in t on the left side 
and, using similar arguments, analyticity in the lower 
half of the z plane goes hand in hand with analyticity 
in t on the right side. We utilize the analytic properties 
of the vacuum waves in space and time as a marker for 
distinguishing the vacuum modes. 



D. Combinations 

The vacuum states of slow-light polaritons are char- 
acterized by analytic wave functions in z. Therefore, to 
describe the vacuum after the formation of the horizon, 
we should construct orthonormal combinations v of the 
non-analytic u and w modes that are analytic on cither 
the upper or the lower half of the complex z plane at 
some arbitrary time to- The solution is 



le 



vr — hsech(nfi)wL — iwR — ie 2tlAjto ^sech(n fi)w E 



R ' 



Vl = ^sech(7r^i)w/j — iw^ — ie~ 



2iwt °isech(7r/iX, 



(54) 



Because the u and w modes are orthonormal with respect 
to the scalar product (pi|), one can easily verify that the 
v modes form an orthonormal set as well. 

The modes ( pi) ) are given on the real z axis and are 
subject to analytic continuation. We prove that the vr 
and modes are analytic on the upper half plane and 
that the vl and v^ modes are analytic on the lower half 
plane. First, we use the definitions (^) and (p0|), and 
write vr and vl in the form 



2cosh(7T/i)vl — e 2v ^vr 
= *£- ie +w >* ((1 + e' 2 ^)u+ 
- e - *" [u+ - ie-*" ((1 ■ 



-27r/j^-* -2iwio 
R 



+2nfj, +* -2iuit 



2cosh(7T/i) vl — e 2ix ^vl 



ie - *" ((1 + e +2nfl )u+ 

: — 2iuto 



Then we show that the combinations 



IE 



(55) 



(56) 



are analytic on the upper half plane and that the corre- 
sponding combinations 



L > U R le 



+*/*«+* 



ie "*u+ , u R 



-2iuto 



(57) 



are analytic on the lower half plane. The analyticity of 
the modes ( [54] ) follows from the analytic properties of 
the combinations (^) and ([57]). Here it is sufficient to 
focus on the vicinity of the origin where the left and right 
modes are connected. As a consequence of Eqs. (|3S|), ( |46| ) 
and (Erf), we get in terms of C = kz, 



e(z) 



«£(0 = «£(-£)• 



(58) 



Consider the analytic properties of Qt 1 -+ 1 / 2 for an arbi- 
trary real \i. We indicate with a suffix ± whether ^v+i/ 2 
should be regarded as analytic on the upper (+) or on 
the lower (— ) half plane, respectively. With this notation 
we get 

= e(c) c v+1/2 + e(-o {-i) l £ +1/2 (_o v+1/2 

= 6(C) (^ l+1 / 2 + 9(— C) e ±l7T ^ A1+1 / 2 ^ (-C)v +1 /2 
= 6(C) C l+1/2 ± 6(-C) (-C) v+1/2 • (59) 

This relation proves the analyticity of the combinations 
(f56l ) and ([57]) and, as a result, the analyticity of the 
modes (j54|). 

However, is the set of modes ( |54"|) unique? In princi- 
ple, we could perform linear canonical transformations 
that convert the modes ( |54|) into a new set, yet such 
transformations are severely restricted by the required 
analyticity in space and time. For example, we can 
construct superpositions of the vr and modes to 
form the new modes vr cos 8 + e^v^ sin 9 and v^ cos 9 — 
e~ ll VRSui9, analogous to the mode transformation of 
a beam splitter f34| . Or we may combine vr with v* L 
in the Bogoliubov transformation vr cosh £ + e ll v* L sinh £ 
and i>^cosh£ + e I7 v£sinh£, analogous to a parametric 
amplifier |35| ]. However, only transformations of this 
type maintain the analyticity on one of the half planes 
of complex z. We are not allowed to combine vr with 
Vft or vr with vl- Yet the possible mode transforma- 
tions are further restricted: a Bogoliubov transformation 
vrcos1i£ + e*^t?2 sinh£ would generate negative frequen- 
cies on the left side of the horizon. We have argued 
that this must not happen. Similar arguments apply 
to all other Bogoliubov transformations. Consequently, 
our modes are uniquely defined up to superpositions, but 
such transformations do not change the vacuum state 
pS4j| . Therefore, the v modes are indeed the vacuum 
modes. 
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E. Radiation 

As a consequence of a slow-light catastrophe, the po- 
lariton field is decomposed into two different sets of 
modes. The v modes contain the polariton vacuum, 
whereas the u and w modes guide the detectable quanta, 

0= [a R w R + d R ±u R ; + a L w L + a L ±u L + H.c. ) dq , 



b R v R + b R ±v R + b L v L + b L1 _vi + H.c. ) dq . (60) 



The a operators are the annihilation operators of the de- 
tector modes and the b operators refer to the vacuum 
modes. H.c. denotes the Hermitian conjugate and the 
q are the mode indices k/ko- Notice that the vacuum 
modes contain both positive and negative frequency com- 
ponents, because the set (|54|) involves complex conjugate 
w modes. Therefore we expect that the corresponding 
operator transformations combine annihilation with cre- 
ation operators. This is the decisive sign of pair creation, 
similar to the production of photon pairs in parametric 
downconversion p2| ]. We represent d R as (w R ,ip), a R ± 
as (u R , (f>) et cetera, use the normalization ( |23| ) and the 
properties of the scalar product (24), and arrive at the 
Bogoliubov transformations [Esf 



d R = isech(7r^)fe i - ib R + ie 2lut ° \sech.{ir ii)b R , 
1 



a R ± = 



b R± + ie^b L± 



(61) 



and at analogous relations for the cll and With- 
out initial probe light injected, the dynamically formed 
slow-light catastrophe will cause spontaneous radiation 
of probe polaritons at a constant rate, because we obtain 
for v-mode vacua 

(a' R (qi) a R (q 2 )) = (d{(q 1 )d L (q 2 )) = n5{q 1 - q 2 ) (62) 

with the average particle number 



1 



(63) 



The radiation energy will be taken from the control beam. 
The initial formation of the horizon is a time-dependent 
process that, therefore, transfers energy to the polariton 
field. Yet, in addition to an initial brief burst of energy, 
the control beam creates a wave catastrophe that would 
force polaritons into a state they cannot occupy. The 
frustrated polariton field reacts in attempting to alter 
the parabolic profile of the control intensity. This process 
takes energy away from the control beam and allows the 
creation of polarition pairs. Pair production continues as 
long as the control beam is not significantly depleted. A 
running wave of control light will produce a steady flow 
of slow-light quanta. 



The slow-light catastrophe generates a maximal parti- 
cle number per mode ( |63| ) of 1/4, which is quite substan- 
tial, considering the fact that bright sunlight with a radi- 
ation temperature of 6 • 10 3 K carries a mere 0.01 photons 
per mode in the optical range of the Planck spectrum. 
However, the photon number (^) is sharply peaked as 
a function of /i and, in any case, our pair-production 
mechanism is restricted to the narrow frequency window 
of EIT |l7j . To maximize the generated quantum radia- 
tion, one should create a situation where /i is near zero 
over an as large as possible spectral range. In terms of 
the experimental parameters, we get in the transparency 
window near ujq, 



2 U 



1/2 



UJ — UJq 



So 



8a 2 u^ 



32tt 2 



(64) 



Pair production occurs on the blue side of the critical 
detuning <5o (for 5 > Sq) where Sq also determines the 
width of the spectrum ( p3|) . The smaller the scale a of 
the group-index profile (|3C|) is, the larger is the criti- 
cal detuning (|64|) and the wider is the spectrum ( |6^ ) 
of particle production. Equations (||) and ( |l6| ) indicate 
that the scale a is small for a steep group- velocity profile 
created by a large control-field gradient. Gravitational 



black holes show a similar behavior |10|. The smaller the 
hole is, the larger is the gravity gradient at the horizon 
and the stronger is the Hawking radiation |lC| gener- 
ated. Returning to our case, the parameter-dependence 
of the particle-number spectrum ( |63| ) underlines the cru- 
cial role of a spatially varying group velocity in creating 
a quantum catastrophe. A mere zero of the group veloc- 
ity would not suffice to produce a measurable radiation. 
The control-field gradient matters. 

Close to the horizon the susceptibility of slow light 
diverges. Yet Nature tends to prevent infinite suscepti- 
bilities: Instead of responding infinitely strongly, opti- 
cal media become absorptive or non-linear. Considering 
gravitational black holes, Nature could prevent the exis- 
tence of true event horizons as well. Here waves are sup- 
posed to shrink in wave lengths beyond the Planck scale 
where the physics is unknown. Hawking radiation seems 
to stem from these extremely shortened waves. This 
trans-Planckian problem |36| was analyzed in theoretical 
toy models of sonic black holes in moving fluids p7| , ^8| . 
Here the inter-atomic distance provides a natural cut-off 
for extremely shortened sound waves. Nevertheless [ [38[ , 
the mere threat of a horizon seems to be sufficient for 
generating Hawking sound. 

Returning to slow light, an EIT medium in linear re- 
sponse is transparent within a narrow spectral window 
around the resonance frequency luq- The spectral width 
of the transparency window is proportional to the inten- 
sity of the control beam u%. Therefore, slow-light waves 
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oscillating at frequencies different from uq are absorbed 
near a node of the control field, unless the medium be- 
comes non-linear. We show in Appendix A that the non- 
linearity of the EIT medium depends on the ratio of the 
probe and control intensities. According to the linear 
optics of slow light discussed in the body of this paper, 
the probe intensity is proportional to the distance z from 
the horizon, as long as kz is small, whereas the control 
intensity grows quadratically. Consequently, at a cer- 
tain distance Zq both fields are comparable in strength. 
Here slow light leaves the regime of linear response. In 
a semi-classical concept of light |39|,^0| quantum fluctu- 
ations are small perturbations of the classical amplitude 
and are subject to a linear theory. If the vacuum state 
of light is classically unstable, photons are created spon- 
taneously. For example, in parametric downconversion 
p2f quantum fluctuations are amplified p^j , generating 
photon pairs. The instability of the linear optics near 
a slow-light horizon may provide the microscopic mech- 
anism for the pair production we have predicted phe- 
nomenologically. However, the non-linear effects of the 
EIT medium are required to dominate at a stage where 
the absorption is still small. We show in Appendix B 
that the probe intensity is proportional to the detuning 
5q given by Eq. (Q) in terms of the characteristic scale 
© of the group index profile (|J). Therefore, the ra- 
tio of the control and probe intensities does not depend 
on a, and, consequently, the non-linearity distance zq is 
independent of the control- field gradient. In order to 
avoid absorption the control should be strong enough at 
zq which requires a steep field gradient. Using the exper- 
imental parameters of Refs. |L3 15 1, the Rabi frequency 
p2| of the control field should grow at least by 10MHz 
per wave length Ao in distance from the horizon, as we 
show in Appendix B. In this case the scale a = 5 x 10 3 Ao 
and the critical detuning 6q = 10 -10 . Appendix B in- 
dicates that millions of photons are generated per sec- 
ond, amounting to a gentle glow perhaps visible with the 
naked eye. 



IV. SUMMARY 

Tuning the control field towards a parabolic intensity 
profile causes a catastrophic situation for slow-light po- 
laritons. In turn, the polariton field sets out to deplete 
the control beam, in an attempt to alter the intensity 
profile that has caused the wave catastrophe in the first 
place, yet in vain. The control beam continuously replen- 
ishes the parabolic intensity profile, driving a stationary 
production of polariton pairs. The two polaritons of each 
pair are created on opposite sides near the horizon, they 
depart at a snail's pace, accelerate gradually and emerge 
as detectable photons. The Hawking radiation of a black 
hole jlO| follows a similar scenario Q. Here the gravita- 
tional collapse 0| has triggered a quantum catastrophe at 



the event horizon, causing pair creation lasting as long as 
the hole possesses gravitational energy |8|JlCt| . One parti- 
cle of each pair falls into the black hole, whereas the other 
escapes into space and appears as thermal radiation flio| . 
In our case, and in contrast to gravitational holes, one 
can explore the other side beyond the horizon and, for 
example, measure the correlations of the generated pho- 
ton pairs. Both cases are triggered by catastrophic events 
with lasting consequences. 

The quantum radiation of a slow-light catastrophe re- 
sembles Hawking radiation but also exhibits some in- 
teresting differences. The emitted spectrum ( |63| ) is not 
Planckian, whereas a black hole of Schwarzschild radius 
r s appears as a black-body radiator with temperature 
hc/(4:TTr s ) juj. The differences between the two spec- 
tra can be traced back to two different classes of wave 
catastrophes. In both cases, waves freeze at an horizon 
in the form £ p with an exponent ifi + 1/2 for slow-light 
media but with an exponent ifi for black holes where 
fi = 2irr s ui/c Note that Unruh's effect (|] of ra- 
diation seen by an accelerated observer is of Hawking- 



class as well 
black holes [2 



and so are most of the proposed artificial 
Remarkably, Schwinger's 
pair production of charged particles in electrostatic fields 
fl4| is accompanied by a subtle wave catastrophe of expo- 
nent ifi — 1/2 Q and leads to a Boltzmannian spectrum 
n = cxp(— 2nfi). All three catastrophes agree in the limit 
of large fi but deviate significantly in the regime of max- 
imal particle production where fi is small. It might be 
interesting to find out whether more than three types of 
quantum catastrophes can occur. 
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APPENDIX A: DARK-STATE DYNAMICS 

In this appendix we consider the microscopic theory 
of the atoms constituting an EIT medium. We derive 
the wave equation (j^) as the linear-response limit of the 
non-linear dynamics of slow light. Assume that the EIT 
medium consists of ua identical atoms per unit volume, 
each one equipped with three levels interacting near- 
resonantly with the probe and control fields respectively, 
see Fig. |[ 
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3 

Probe 




Control 



FIG. 3. Three-level atom in a regime of Electromagnet- 
ically-Induced Transparency. The control beam couples the 
levels 2 and 3, which influences strongly the optical properties 
of the atom for a weaker probe beam tuned to the transition 
l<->3. 



We treat the light as a classical electromagnetic field. 
An atom is characterized by the energy-level differences 
Huji2 and hu>23 with w 12 + w 2 3 = W13 = ^o- Typically, the 
transition frequencies W13 and U23 are in the optical range 
of the spectrum or in the near infrared (10 15 Hz), whereas 
the frequency u>±2 is much lower (10 9 Hz). The atom is 
subject to fast relaxation mechanisms (10 6 Hz) that trans- 
port atomic excitations from the | 3 ) state down to | 1 ) 
and from | 3 ) to | 2 ) , mainly caused by spontaneous emis- 
sion. Hardly any excitations move from | 2) to | 1), be- 
cause the spontaneous emission rate is proportional to 
the cube of the frequency pj| . Here the relaxation may 
be dominated by other processes, for instance by spin- 
exchanging collisions. Without relaxation the dynamics 
of the atom is governed by the Hamiltonian 



H = 






1 K + ) 





hu\2 

h K23 E, 



1 «■ z?( - ) 



(+) 



13 



(65) 



The Hamiltonian represents the atomic level structure 
and describes the dipole interaction with light, consid- 
ering here only the positive/negative frequency compo- 
nents Ep and Ec that match approximately the level 
structure. The E p and E c fields are the probe and con- 
trol light respectively and are given in SI units. We de- 
scribe relaxation phcnomcnologically by the transition 
processes 



.4, 



(66) 



occurring at the rates 71 and 72, typically a few 10 6 Hz. 
The density matrix of the atom, p, evolves according to 
the master equation Efl] 



dp * r~ rVi 



(A\A lP - 2A t pAj + pAjAt) . (67) 

It is advantageous to represent the light fields in terms 
of Rabi frequencies 



EW 



£2 r, 



-iujQt 



Kl3 



E (+) 

h p 



(68) 



defined here with respect to the atomic transition fre- 
quencies oj c = W23 and ujq = W13. In the absence of relax- 
ation, an atom would oscillate between the ground and 
the excited state with frequency SI (Rabi flopping p2|). 
On the other hand, relaxation leads to a stationary state 
where the atomic dipoles follow the fields. 

Assume that the control beam is in exact resonance lu c 
and that the probe light is monochromatic with a small 
detuning lo — luq. Furthermore, the Rabi frequency of the 
control beam shall dominate all relevant time scales, 



|fi c | > I O p | , 71, 72,|w-w c | 



(69) 



In this limit the stationary state of the atomic evolution 
( |67| ) turns out to approach the pure state 



P=\ipo)(ipo 
I ) = U Q N [ I 



£2p 



2 ^ ~ g°) n 13' 



(70) 



that is called a dark state p7|| . Here we have separated 
the rapid oscillations of the atom at the optical transition 
frequencies from the slower atomic dynamics, 



U 



1 

e~ luJl2t 







-iujQt 



(71) 



Suppose that a dominant and monochromatic control 
beam has, after relaxation, prepared the atom in the sta- 
tionary state (|7(]). How will the atom evolve when the 
control and probe strengths vary? First we show that the 
atom remains in a pure state, as long as the | 3 ) compo- 
nent is small, 



P 33 = (3|/5|3)«l 



(72) 



Consider the statistical purity tr{p 2 }. A quantum system 
is in a pure state if and only if the purity is unity . We 
apply the master equation ( |67j ) and see that the purity 
does not change significantly, 

dtr{p 2 } = 2tr{pdp} 

= 4 [71(1 - p n ) + 72(1 - P22)] P33*, (73) 

once the atom has occupied a pure state with sparsely 
populated level | 3 ) . Consequently, we can describe the 
state of the atom by a vector | ip ) . 
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Suppose that the control and the probe strengths vary. 
How does a dark state follow the light? In the case ( |7^ ) 
the state vector is dominated by its components in the 
subspace spanned by the two lower levels 1 1 ) and | 2 ) . If 
we find a vector | ip ) that describes correctly the dynam- 
ics ( |67j ) in this subspace, the third component ( 3 | ip ) 
must be correct as well, to leading order in p^. The 
lower ranks enslave the top level. Since the relaxation 
processes (^6|) do not operate within the lower subspace, 
we can ignore dissipation entirely, to find the dominant 
state of the atom. We write down the state vector 



\rl>) = U N[\l) 



n r 1 ' 



with the abbreviations 



a 



No=[l 



jo 



%\ 2 



-1/2 



N = N exp 



n p \ 2 de 



(75) 



In a stationary regime under the condition ( p9| ) the vector 
(f74|) agrees with the dark state (f7(i|). We see from the 
properties 

d t N = -NN*-£df£, d t N-£ = NN*d t -£ (76) 
that | ip ) satisfies the differential equation 



ihdtlip) = H\ip) +ihd t (3\ip)\3] 



(77) 



Consequently, the vector ( f74| ) describes correctly the dy- 
namics of the atom in the lower-level subspace. There- 
fore, the atom remains in the dark state (]74|), as long 
as the atom's evolution never leads to an overpopulation 
at the top level | 3 ) . The initial relaxation-dominated 
regime has prepared the dark state, but later the atom 
follows dynamically without relaxation p9| . 

In responding to the light fields, the evolving atoms 
constitute a macroscopic dipole density called the matter 
polarization P A . Consider a one-dimensional model for 
light propagation. The matter polarization influences the 
probe light according to the wave equation 



(78) 



Each atom generates a dipole moment 
of 5Ki3tr{p| 1)(3 | }, oscillating at positive frequencies, 
that contributes to the total dipole density. Therefore, a 
medium with n A atoms per volume generates a matter 
polarization with the positive-frequency component 



P 



(+) _ nA 



^13 (3 1 -0>( ^ 1 1 ; 



n A ■ 



K 13 JV 



n \n r \ 2 



x I idt — ujq — i 



.(dt\Sl c \) 



(79) 



where 9 C = argf2 c . Assume, for simplicity, that Q c is real. 
Otherwise we can easily incorporate the phase 9 C of the 
control field in the phase of the electric field without af- 
fecting the wave equation ( |78| ) , as long as 9 C varies slowly 
compared with the optical frequency loq. We define 



a = n A 



ojq n A Tiuq 



2e ?i I O c | 2 2 k\ 2 s \E c \ 2 



(80) 



which, as we will see shortly, is the group index (Q). We 
get 



e 1 5 t 2 P« w -N*a 2w (id t - wo - 



Et 



51) 



and approximate 

2io (id t - lo ) 4+) « {idt + Uo)(idt - lo ) 4+) 

= -(a?+w 2 )M+). (82) 



In this way we obtain from the general wave equation 
(f78| ) an equation that is valid for both the positive and 
the negative frequency component of the probe light, 



[8 2 t - c 2 dl + 7V 4 (dta d t + aul)]E p = Q. 



(83) 



The dark-state dynamics may lead to a non-linear sat- 
uration of the medium, described by the Nq factor in 
the wave equation (|S3"|). The non-linearity is relevant 
when the Rabi frequencies \Q P \ and |f2 c | are compara- 
ble. When the probe is significantly weaker than the 
control light, the medium responds linearly, 



(d t (l + a)d t 



auj n 



E n = 0. 



(84) 



We have derived the wave equation (|^). The group in- 
dex (|80| ) determines the group velocity v g = c/(l + a). 
Remarkably, the lower the intensity of the control beam 
is, the slower the probe light becomes. Taken to the ex- 
treme, light freezes when the control light is switched 
off — a paradoxical behavior that is only possible in a 
dynamical regime ps|| : a control beam of moderate in- 
tensity first captures the probe light, slowing it down, 
and then, by ramping down the control intensity, freezes 
the probe pulse. Equally paradoxically, the non-linearity 
of the EIT medium is stronger the weaker the control 
beam is. We show in Appendix B that the unusual non- 
linear optics in an EIT medium matters in a slow-light 
catastrophe. 
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APPENDIX B: ESTIMATIONS 

In this appendix we estimate the effect of a slow-light 
catastrophe, using the experimental parameters of Refs. 
jl3],[l5|]. We calculate the photon flux in the far field 
and estimate the intensity near the horizon. In order 
to detect experimentally the quantum radiation of the 
catastrophe, the flux should be sufficiently strong. Close 
to the horizon, even slightly detuned slow light leaves 
the absorption-less transparency window of EIT, unless 
the medium becomes non-linear. The intensity near the 
horizon determines whether the non-linearity or the ab- 
sorption dominates. 

First we generalize our one-dimensional model of the 
slow-light catastrophe to the three dimensions of space in 
Cartesian coordinates x = (x, y, z). The spatial profile of 
the group index shall be uniform in x and y and parabolic 
in z with the scale a. The propagation of slow light is 
governed by the wave equation 



d 2 - c 2 v 2 



ip = 0. 



(85) 



We find the stationary solutions and normalize them ac- 
cording to the scalar product 

/2 
(ip*dttp2 - tpid t ip*^(l + % I )d 3 x. (86) 

We obtain the set of modes (pOf) with 



Q(z) 



A V2c 
1 
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w-j 



wf (x) =w%(-x), 



(87) 



and Ao = 2ir/ko and ko — loq/c. Armed with the three- 
dimensional modes, we turn to calculating the energy 
flux. The Poynting vector of light is the time-averaged 
expectation value of the normally ordered Poynting op- 
erator Q, 

P = ~L [ +T (-- -hc 2 (d^)(V£) :) dt 

I — >OC 11 J _rp 

(: W1WW2 ai&l + w\S/w2 a\a,2 :) d 3 qi d 3 q2 



-he 2 



-he 2 



(wVw* + lii'Vw) n d 3 q 



he J 2uj 2 q\w\ 2 nd 3 q . 



(88) 



In the last step we have utilized that the w modes ap- 
proach plane waves @ far away from the horizon. In 
view of the narrow bandwidth of EIT we can replace uj 



by the resonance frequency loq. Consider a radiating sur- 
face with area A observed from the distance r under the 
angle <&. The range of wave vectors contributing to the 
flux (§8|) is restricted to lie within the solid angle of the 
surface, Acos'd/r 2 . In the line of sight we thus get the 
Poynting-vector component 

j4cos?9 f 2 - . 

F = hujQ j — 2cujq / | to I ndq 



A cos i? 
ftw ^ — uj q 



1 - e- 27rtl 



2tt 



dq, 



(89) 



having applied the asymptotics (|5l]). We employ \i as the 
integration variable, with fi(d/.i/dq) = a 2 ko, and obtain 
in terms of the critical detuning (64) the photon flux 
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(90) 



The flux integrated over the two half spheres around the 
radiating surface gives the total photon-production rate 



W 2 p r 2 j{ 
N = Air - — sin $ d$ = —loqSq 2Trrjo . 
Jo 



hujo 



\2 
A 



(91) 



For light in the optical spectral range, 27rwo Vo is about 
4 x 10 14 Hz. Assuming a critical detuning So of 10~ 10 
a surface of 10 2 Ao could generate millions of photons 
per second. Usually a photodetector captures only a 
small solid angle of the radiation emitted from a local- 
ized source, and a detector is not perfectly efficient in 
counting all photons. Yet the radiation of the slow-light 
catastrophe seems to be strong enough to be detectable. 

Let us estimate the strength of the control field needed 
to generate the flux we have calculated. We calibrate the 
field strength in terms of a Rabi frequency (|8|). Ac- 
cording to Eq. (U), a critical detuning Sq of 10~ 10 cor- 
responds to a length scale a«5x 10 3 Ao of the group- 
index profile (^0|). In the experiment Jll| a group index 
of 10 7 is generated by a control field with Rabi frequency 
Q c = 2.57 x 27rMHz w 16MHz. The group index is in- 
versely proportional to Jl 2 , which leads to a Rabi fre- 
quency fl c of about 5 x 10 4 v / aMHz. For the profile ( |30| ) 
of our slow- light catastrophe we find that Q c should grow 
linearly by 10MHz per wave-length distance away from 
the horizon, 



n c = 10 — MHz . 
A 



(92) 
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Probably 10 to 100 wave lengths are sufficient to establish 
the wave catastrophe. This would take a continuous- wave 
control field with a maximal Rabi frequency of 10 2 to 
10 3 MHz. 

The ratio £l p /Q c of the Rabi frequencies j68| ) deter- 
mines whether non-linear effects dominate near the hori- 
zon. We calculate £l p using the results of the linear the- 
ory. Strictly speaking, Rabi frequencies refer to classical 
fields. Here we regard the time-averaged and normally 
ordered expectation value of the quantum intensity E p 
as being proportional to il p . Comparing Eq. (1) of Ref. 
|L5j with our Eq. ( |S0| ) we find a relation between the Rabi 
frequency and the intensity 



3x 10 3 (^)m 3 Hz 



(93) 



for the atomic transition employed in the experiments 
jl3|,[l5| . We express E p in terms of the field (p in units of 
the vacuum noise, and get 



hajQ ' 



2uj J \w\ 2 nd 3 q. 



(94) 



Close to the horizon the w waves obey the asymptotics 
, l9 koz 1 



w 



2c\l 1 



(95) 



as we find from the definition ( p7| ) and the behavior (pq) 
of the Bessel functions, utilizing Eq. 1.2.(6) of Ref. |33] ] 
for the gamma function x\ = T(x + 1). We obtain 



I |2- ^ 
\w\ n 



a 2 k% 



dfi 



327r(l-ln2) l 5 



*8 



and, consequently, 



n p ~ 8.6 W— MHz. 



(96) 



(97) 



Therefore, the Rabi frequencies tt p and £l c are compara- 
ble at about half a wave-length distance away from the 
horizon. Here a critical detuning of 10~ 10 lies still within 
the transparency window generated by a control field of 
5MHz Rabi frequency. The EIT medium becomes non- 
linear. 
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